Overlaying optical lattices for simulation of complex frustrated antiferromagnets 
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We present design techniques of special optical lattices that allow quantum simulation of spin 
frustration in two-dimensional systems. By carefully overlaying optical lattices with different periods 
and orientations, we are able to adjust the ratio between the nearest-neighbor and next-nearest- 
neighbor interaction strengths in a square spin lattice and realize frustration effects. We show 
that only laser beams of a single frequency is required, and the parameter space reachable in our 
design is broad enough to study the important phases in the J1-J2 frustrated Heisenberg model and 
checkerboard antiferromagnet model. By using the polarization spectroscopy for detection, distinct 
quantum phases and quantum phase transition points can be characterized straightforwardly. Our 
design thus offers a suitable setup for simulation of frustrated spin systems. 
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INTRODUCTION 



Geometrically frustrated systems with strong correla- 
tion have attracted much attention due to the highly 
nontrivial interplay between frustration and correla- 
tion 0, || • In such systems, not all pairs of spins can 
simultaneously assume their lowest energy configuration, 
and the ground state usually has a large number of de- 
generacies. In particular, frustration in quantum antifer- 
romagnets may cause certain types of magnetically disor- 
dered quantum phases, including the resonating valence 
bond spin liquid state Q and the valence bond crystal 
state As an exotic state of quantum matter, the spin 
liquid in quantum frustrated antiferromagnets is a topic 
of great interest in condensed matter physics. In such 
states, the local moments do not form ordering down to 
the very low temperatures despite of strong antiferromag- 
netic couplings. It has been proposed that these exotic 
states are closely related to the mechanism of high-T c 
superconductivity For instance, the quantum spin 
liquid state could become unconventional superconduct- 
ing state when the charge carriers are introduced. 

Despite intense studies over the past several decades, 
the nature of spin liquid in the strongly frustrated regime 
remains poorly understood. A typical system is the frus- 
trated J1-J2 Heisenberg model on a square lattice. Be- 
sides the nearest neighboring (nn) spin interactions lead- 
ing to Neel state, the J1-J2 Heisenberg model contains 
extra next nearest neighboring (nnn) spin interactions. 
Frustration may lead to the destruction of long-range or- 
der like Neel order, and quantum disordered phases may 
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show up in the strongly frustrated regime with apprecia- 
ble value of Jij J\. Many candidates for the ground states 
have been suggested, such as a columnar spin dimer- 
ized state, a plaquette resonating valence bond phase, 
and a columnar spin dimerized state with plaquette-type 
modulation. Because of the difficulties in theoretical 
and experimental studies of the frustrated spin liquid 
physics, a quantum simulator for such systems is highly 
valuabledi]. 

Recent advance in the manipulations of ultracold 
atoms in optical lattice has opened new possibilities for 
realizing simulation of many significant strongly corre- 
lated quantum models @, Q- Much efforts have been 
devoted to the simulation of quantum magnetism in ul- 
tracold lattice gases. Recent experiments have shown 
the evidences of superexchange interactions between two 
neighboring sites [8|. Very recently, the realization of 
large scale quantum simulator of frustrated magnetism 
in triangular optical lattices has been achieved Q. As 
for the simulation of the geometrically frustrated antifer- 
romagnets, an effective manipulating method is highly 
demanded to reach the strongly frustrated regime of the 
J1-J2 Heisenberg model where the exotic quantum spin 
liquid states may appear around J\j J% ~ 0.5. However, 
it is technically difficult to achieve the interesting physi- 
cal regime. As we know, the amplitude of intersite tun- 
neling rate decreases exponentially with their distance 
and thus the nnn tunneling rate is much smaller than 
the nn tunneling rate in a normal square lattice. Ac- 
cording to the mechanism of superexchange interactions, 
the value of J2/J1 is much more suppressed. 

Another outstanding difficulty in exploring magnetic 
systems in cold atoms is the detection of quantum many- 
body states. Various detection methods have been pro- 
posed during the past years. The noise correlations tech- 
niques El may detect density-density correlations, 
which corresponds to the spin-spin correlations. Anal- 
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ogous to the neutron diffraction in condensed matter 
systems, Bragg scattering spectroscopy [l2[ may probe 
the static spin structure factor of the Neel state. Single- 
lattice detectors [l3| may measure magnetic ordering, but 
is still a kind of destructive technique. Recently, the po- 
larization spectroscopy technology [lH was proposed to 
detect magnetic correlations including ferromagnetic and 
antiferromagnetic but also more exotic spin ordering by 
tuning the parameters of polarized light, which makes it 
a valuable tool for quantum states measurement. 

Engineering of optical lattices has proven very versa- 
tile as shown in previous studies [HI, [l6| where suitable 
laser beam configurations have been employed to pro- 
duce sophisticated optical lattices. In this work, we push 
the idea further by overlaying optical lattices with dif- 
ferent periods and orientations to realize spin frustra- 
tion effects in two-dimensional square lattices. Filling 
such a carefully designed optical lattice with ultracold 
fermionic atoms, we can then manipulate spin-spin ex- 
change interaction between two next-nearest neighboring 
atoms in optical lattice by adjusting the optical poten- 
tials. We calculate the tunneling rate in the square lattice 
and find that the nnn tunneling rate could be compara- 
ble to the nn tunneling rate. It turns out that the J\-J<i 
Heisenberg model in the strongly frustrated regime could 
be realized in which the exotic spin liquid states may 
exist. Moreover, we can simulate other geometrically 
frustrated lattice like the checkerboard antiferromagnet 
model. As for the detection of quantum many-body 
states, we utilize the polarization spectroscopy technol- 
ogy [3 to study the spatially resolved spin-spin correla- 
tion functions. For J\-Ji Heisenberg model, we demon- 
strate that this promising method can be used to de- 
termine the quantum phase transition points straight- 
forwardly and characterize distinct quantum phases. In 
particular, one-to-one correspondence between the ex- 
treme points of noise signal and quantum phase tran- 
sition points has been achieved. 

The rest of the paper is organized as follows. In Sec. 
II, we implement the two-dimensional frustrated optical 
lattice using a special optical engineering technique. In 
Sec. Ill, simulations of the J\- J2 Heisenberg model and 
the checkerboard antiferromagnet model are introduced. 
In Sec. IV, we employ the polarization spectroscopy tech- 
nology to detect the quantum many-body states. Finally 
Sec. V gives a summary. 



II. THE TWO-DIMENSIONAL FRUSTRATED 
OPTICAL LATTICE 

The Heisenberg spin interactions between neighboring 
ultracold atoms in an optical lattice arise due to the in- 
tersite virtual tunneling of atoms [H, Q3] • Consequently, 
the spin-spin exchange interaction depends on the inter- 
site tunneling rate t as well as the on-site Coulomb inter- 
action U via the superexchange mechanism, J ~ t 2 /U. 
To implement a J\-J% Heisenberg model in the strongly 



frustrated regime, the crucial procedure is to design a 
special two-dimensional optical lattice in which the nnn 
tunneling rate t% can be comparable to or even larger 
than the nn tunneling rate t\. It is hard to achieve by 
creating ordinary optical lattices using laser beams of the 
same wave number k as is required to trap atoms in a cer- 
tain state. However, by properly designing several optical 
lattices with laser beams traveling in different directions 
and superimposing them elaborately, we can tune the ef- 
fective potential barriers along different directions and 
realize such a frustrated optical lattice using laser beams 
with the same wave number k. 

We start by creating a normal two-dimensional square 
lattice. First, we use a strong standing-wave field in the 
z direction to completely suppress the tunneling of atoms 
along the z direction. Then, as shown in Fig. 1(a), we 
shine a pair of blue-detuned interfering traveling laser 
beams of wave number k in the z-x plane, each with an 
angle of 7r/4 with respect to the x-y plane. These two 
laser beams co-propagate in the z direction and counter- 
propagate in the x direction. Together with the strong 
trapping field in the z direction, it creates a periodic 
potential along the x direction in the x-y plane. Sim- 
ilarly, as shown in Fig. 1(a), we shine one more pair 
of blue-detuned interfering laser beams in the y-z plane, 
each with an angle of 7r/4 with respect to the x-y plane. 
This creates a periodic potential along the y direction in 
the x-y plane. The net effect of these trapping lasers is 
the formation of the following two-dimensional spatially 
varying potential profile: 

V x {x,y) =Vi[sin 2 (^=x+^) + sin 2 (^=y + ^)}, (1) 

where k is the wave vector for the laser beams and V± 
is proportional to the dynamical atomic polarizability 
times the laser intensity. Its projection into the x-y plane 
is k/V2. This potential profile then produces a square 
optical lattice with a period a\ = \j \[2 where the wave- 
length A = 2ir/k. We load the optical lattice such that 
there is one atom trapped at each minimum of the optical 
potential. In such an optical lattice, since the tunneling 
rate is sensitive to the potential barrier height and lat- 
tice site separation, it can be easily seen from Fig. 1(b) 
that the nnn tunneling rate along the diagonal direction 
is much smaller than the nn tunneling along the x, y di- 
rections. This is because the potential barrier along the 
diagonal direction is higher and the distance between the 
sites in that direction is also larger. 

In order to make the nnn tunneling rate ti comparable 
to the nn tunneling rate t\ and thus realize an appreciable 
J2/J1, we must adjust the potential barrier heights along 
both the diagonal and x, y directions. For this purpose, 
we add another square lattice along the diagonal direc- 
tion in the x-y plane with a lattice period 0,2 = A/2. As 
shown in Fig. 1(c), this second lattice can be created by 
applying two pairs of standing-wave laser beams in the 
x-y plane along the diagonal directions and of the same 
wave number fc as those used earlier to engineer the po- 
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tential profile in Eq. 
form 



([T]). Its lattice potential has the 



k 



V 2 (x,y) =V 2 [sin {-j={x + y) 



sin 2 {-j={x-y))\. (2) 



To avoid undesired interferences between laser beams, 
one can randomize the relative orientation of the polar- 
ization of laser beams generating optical lattice poten- 
tials of Vi and V 2 [15]. Thus, the total effective two- 
dimensional optical potential is a sum of Vi(x,y) and 
V 2 (x,y). Its profile is plotted in Fig. 1(d) for V% = 10E R 
and V 2 = 15E R where E R = h 2 k 2 /2m is the atomic re- 
coil energy. The corresponding lattice configuration is 
depicted in Fig. 1(e). It is clear that due to the sec- 
ond lattice we engineered, the effective potential barriers 
along x, y directions are much enhanced and the tunnel- 
ing between nearest neighbors is suppressed. In contrast, 
the effective potential barrier along the diagonal direction 
is lowered and the tunneling between nnn is enhanced. 
Thus, we can tune the ratio t 2 /ti by adjusting the value 
ofVa/Vi- 

To find out more practical parameter region of t 2 /ti in 
a physical system, we numerically calculate these tunnel- 
ing rates for an optical lattice loaded with Sodium atoms. 
The tunneling rate between site i and site j is shown as, 



U 



(WilHolWj, 



(3) 



where \Wi), \Wj) are the Wannier functions [lg, [19( at 
sites i and j, and H = — V 2 /2 + Vi(x, y) + V 2 (x 7 y) de- 
notes the single-particle Hamiltonian of the system. Here 
we consider the lowest band only. For the Sodium atoms, 
Er/H = 2tt x 32 kHz for a blue detuned A = 514 nm. We 
may choose V\ — 10E R and increase V 2 from to 2QE R 
and then plot t\ and t 2 in Fig. 1(f). It is obvious that 
the value of t 2 (ii) increases (decreases) as V 2 increases. 
Notably, t 2 becomes larger than t\ when V 2 > 15E R . 
Therefore, by using a sufficiently large V 2 , we can simu- 
late the J\- J 2 Heisenberg model in the strongly frustrated 
regime. It is worth to mention that we have shown the- 
oretically the way to simulate an one-dimensional J\-J 2 
model in concatenated microcavities [20]. However, it is 
quite difficult to realize a two-dimensional J\-J 2 model in 
the strongly frustrated regime by employing this method. 

Moreover, we can realize a checkerboard lattice by 
adding a third square lattice with a lattice period a 3 = X 
along the diagonal direction in the x-y plane. Since the 
period of this third lattice is twice of the second lattice 
and both are oriented in the same direction, the diagonal 
tunnelings in half of the lattice are strongly suppressed 
due to the presence of potential maxima, and a checker- 
board lattice can be constructed. To generate the third 
lattice, we may adopt a technique similar to that used for 
the first simple lattice. In Fig. 2(a), m is one diagonal 
direction in the x-y plane and m! is the other diagonal 
direction perpendicular to m. We shine a pair of blue- 
detuned interfering traveling laser beams of wave number 
k in the m-z plane, each with an angle 7r/3 with respect 
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FIG. 1: (Color online) (a) Schematic diagram for generating 
the interference field Vi(x,y). (b) Normal square lattice pat- 
tern induced by the field Vi(x, y) in unit of E^. (c) Schematic 
diagram for generating the interference field Vz(x, y). (d) The 
square lattice pattern induced by the field V\ (x, y)+V2(x, y) in 
unit of En. (e) Schematic diagram of the square lattice with 
diagonal tunneling, (f) The nearest-neighboring tunneling ti 
(black solid line) and the next-nearest-neighboring tunneling 
t2 (red dashed line) as a function of V2 in unit of Er. 



to the x-y plane. Similarly, another pair of blue-detuned 
traveling laser beams of the same wave number are shot 
in the m'-z plane, each with an angle 7r/3 with respect to 
the x-y plane. Obviously, these two pairs of laser beams 
will form a two-dimensional lattice along the diagonal di- 
rections in the x-y plane and its period is determined by 
the projected wave number in the x-y plane k^ = k/2. 
Specifically, the spatially varying potential profile for the 
third lattice is 



V, 



\[2 \[2 
,ix,y) =V 3 [sin 2 (—k(x + y)) + sin 2 { — k(x - y))] 



(4) 



For amplitudes V x = 10E R , V 2 = 10E R and V 3 = 5E R , 
the superposed potential V(x,y) — Vi(x,y) + V 2 (x,y) + 
V 3 (x,y) is plotted in Fig. 2(b). We can see that there 
are two sets of plaquettes in which the potential barrier 
heights along the diagonal direction are raised and low- 
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FIG. 2: (Color online) (a) Schematic diagram for generating 
the interference field V-i(x,y). (b) The optical lattice induced 
by interference field V(x,y) in unit of Er. (c) Schematic 
diagram of a checkerboard lattice. 



ered respectively. Consequently, the diagonal tunnelings 
in these two sets of plaquettes are either suppressed or 
enhanced, providing us the checkerboard lattice. The 
configuration of the checkerboard lattice is shown in Fig. 
2(c). When we fix the value of V\ and V3 and adjust 
the value of V2, the rate of ^Ai will be changed accord- 
ingly. When £2 is comparable to £1, the checkerboard 
J\-J% model may reach the strongly frustrated regime 
J 2 /Ji~l. 



III. THE FRUSTRATED Ji-J 2 HEISENBERG 
MODELS 



We can implement the J\-Ji Heisenberg model by us- 
ing a collection of ultracold fermionic atoms confined in 
the special square optical lattice at integer filling. In the 
real experiments, ultracold fermionic atoms have been 
used to realize the Mott insulating state [U [22]. By 
further lowering the temperatures, two easily accessible 
hyperfine states in such system can be used to realize 
magnetic orders 0, [24j . Here we consider that the sys- 
tem enters into the Mott insulating state, and then load 



the ultracold fermionic atoms onto the optical lattice so 
that there is only one atom per lattice site. The two hy- 
perfine states represent the effective spin states \a) =\\) 
and 1 1). We assume that the optical potential has no 
state dependence and the atoms are confined to the low- 
est Bloch band. The model Hamiltonian can be written 
as. 



H = - tx ( C *W + h.c.) 



(5) 



where c; CT is the Fermion annihilation operator for the 
atom on site i, ni a = c\ a Ci„, t\ and £2 are respectively 
the tunneling matrix elements between two nn sites 
and two nnn sites and U denotes the on-site re- 

pulsive interaction between the atoms. The interaction 
U can be adjusted by Feshbach resonances method. In 
the Mott insulator regime, we have £1, £2 <C U. By treat- 
ing the tunneling as weak perturbation and using the 
Schriffcr- Wolff transformation, it can be shown that Eq. 
(O is equivalent to the following J\-Ji Heisenberg model 
Hamiltonian, 



H 



{i,3) 



S, + J: 



(6) 



where Ji = 4£?/{7 (i=l,2), Ji and J2 correspond to the 
nn and nnn exchange couplings, respectively, (see Fig.lc) 
Here we ignore the higher order corrections for the coef- 
ficient Jj, in which the dominant item is proportional to 
tf/U 3 [25[ and can be suppressed greatly as U/U grows. 

Due to the competition between J\ and J2, the frus- 
trated J1-J2 antiferromagnet has certain limiting cases. 
At J2/J1 <C 1, the model is reduced to the Heisenberg 
model on a square lattice. At J2/J1 ~ 1, the system has 
a collinear antiferromagnetic quasi-long range order. In 
the strongly frustrated regime J2/J1 ~ 0.5, there are sev- 
eral candidates of possible spin liquid states, but there 
is still lack of general consensus. According to previ- 
ous analytical and numerical studies [26|-|28[ , there might 
be three quantum critical points around J2/J1 ~ 0.38, 
0.5, and 0.6, respectively. In particular, the simple Neel 
state corresponds to the parameter region Jij J\ < 0.38, 
the simple columnar dimerized spin liquid shows up for 
0.38 < J2/J1 < 0.5, the columnar dimerized spin liquid 
with plaquette type modulation may exist in the param- 
eter region 0.5 < J%l J\ < 0.6, and the collinear state 
may appear when J2/J1 > 0.6. As illustrated in Fig. Id, 
we observe that £2/£i or J2/J1 increases progressively as 
we increases Va- It is clear that the strongly frustrated 
regime of 3\-Ji model can be reached. This model also 
exhibits rich dynamical behaviors which makes it an ideal 
testing ground for the study of quantum phase transi- 
tions. Once the quantum many-body states could be 
faithfully probed, our proposed quantum simulator may 
resolve many outstanding issues. 
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Another interesting frustrated two-dimensional quan- 
tum spin model under our consideration is the checker- 
board antiferromagnet. Similar to the case of J\-J<z 
model, we can simulate the checkerboard antiferromag- 
net model by using ultracold fermionic atoms confined 
in the checkerboard optical lattice. The model Hamil- 
tonian contains both the nn couplings J\ and the nnn 
diagonal link J 2 . In the limit of J%/J\ <C 1, its ground 
state has the Neel long-range order, whereas it corre- 
sponds decoupled Heisenberg chains when J%l J\ 1. 
According to previous studies [HI, [30l | , there may exist 
three distinct quantum phases: the Neel order shows up 
for J%l J\ < 0.8, a valence bond crystal in singlet plaque- 
ttes may appear when 0.8 < Jij J\ < 1.1, and decoupled 
Heisenberg spin chains corresponds to the parameter re- 
gion J2/J1 > 1.1. The checkerboard antiferromagnet 
model in the strongly frustrated regime J2/J1 ~ 1 can 
be realized experimentally with ultracold atoms in opti- 
cal lattices. 



IV. DETECTING THE QUANTUM 
MANY-BODY STATES 

Besides the preparation of the two-dimensional frus- 
trated spin models in optical lattices, another crucial 
procedure is to detect faithfully the quantum many- 
body states in these strongly correlated systems. Among 
many different experimental methods, one of the signifi- 
cant progresses in this field is the quantum polarization 
spectroscopy of the atomic ensemble, which is a type of 
quantum non-demolition detection [14j . After a polar- 
ized probe laser traveling through the atomic ensemble, 
the spin fluctuation and correlation of the atomic sys- 
tem will be mapped on those of light, which can then 
be efficiently measured using homodyne detection. This 
method can be used to probe the spatially resolved spin- 
spin correlation functions and then may distinguish dif- 
ferent magnetic quantum phases. Recently, this method 
has been applied to study the phase diagram of spin-1 
bilinear-biquadratic model and appears to be quite suc- 
cessful 31]. In the following, we will employ the quantum 
polarization spectroscopy method to detect the quantum 
many-body states of a small-scale J\-J-x model and gain 
some important insights on distinct quantum phases and 
quantum phase transition points. 

Next we consider a probe standing wave polarized in 
the x direction characterized by the Stokes operators s\, 
S2, S3. We define a macroscopic Stokes operator S\ = 
J Sidt and introduce the canonical quadrature operators 
X = y/2/ NpS 2 , where Np is the number of photons. 
The variance of X of the output light contains the atomic 
spin-spin correlations: 

((AA ouf ) 2 > = \ + ^~ A Y. c ^(S k z S l z ) (S k z )(S% 

(7) 
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FIG. 3: The noise signal A as a function of J2 / Ji of the J1-J2 
frustrated spin model. 



where k is the coupling constant, F is the total angular 
momentum of a atom, Na is the total number of atoms, 
denotes the z component spin operator of the fcth 
atom. On describes the atom-light coupling and is given 
by 

c n (k p ,b) = 2 J dr cos 2 (k p {r - b))\w(r - r n )\ 2 , (8) 

where k p is the wavevector of the probe laser, b is its shift 
with the optical lattice and w(r — r n ) is the Wannier 
wavefunction of the atom confined at site r„. Here we 
choose F =1/2 and K = 1. The probe laser propagates 
along the diagonal direction with k p = k/A and 6 = 0, 
so the modification of the atom-light coupling c n has the 
structure 

/2 1 1\ 
1012 

12 1' { ' 

\1 2 1 0/ 

The noise signal is defined as A = ((AA OMt ) 2 ) — 1/2, 
which is a linear combination of various spatially varying 
correlation functions. By tuning the parameter like J%j J\ 
of the Hamiltonian, the ground state wavefunction may 
evolve accordingly. When the tuning parameter crosses 
the quantum phase transition point, a sudden change of 
the symmetry of ground state wavefunction may occur. 
The quantum phase transition may lead to the possible 
discontinuity of correlation function itself or its first or- 
der derivative. We naturally expect that the features 
of noise signal as a function of parameter should have 
certain connections with distinct quantum phases and 
quantum phase transition points. Then we numerically 
investigate the evolution of A as a function of parameter 
J2/J1 of the frustrated spin models and expect to reveal 
the connections between noise signal and quantum phase 
transitions. 

We numerically calculate the ground state many-body 
wave function of the spin model by exact diagonaliza- 
tion method for a small cluster 32]. Then we compute 
the spatially varying spin-spin correlation functions and 
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FIG. 4: (a) A as a function of J2/J1 of the checkerboard 
antiferromagnet model; (b) Its first derivative versus J2 / Ji . 



finally we obtain the value of the noise signal. For the 
J1-J2 frustrated model, the noise signal A is plotted as 
a function of J%l J\ in Fig. 3. It can be clearly seen 
that there are two local extreme points at J2/J1 ~ 0.38 
and 0.62, respectively. By analyzing the first derivative 
of noise signal A as a function of J2/J1, we may lo- 
cate another extremum point at J2/J1 = 0.5. According 
to our natural expectation, we may connect these three 
extreme points to three possible quantum phase tran- 
sition points. Very interestingly, these extreme points 
coincide quite precisely with the known quantum phase 
transition points (in an infinite system) respectively at 
{J 2 /J l ) c = 0.38, 0.5 and 0.6 (26H28|. This one-to-one cor- 
respondence between the extreme points of noise signal 
and quantum phase transition points is remarkable. Ac- 
cording to previous studies, two points at ( J2/ J\) c = 0.38 
and 0.6 are connected to the second order quantum phase 
transition while the point at {^/Jijc — 0.5 may cor- 
respond to first order quantum phase transition. That 
may explain why the two extreme points show up for 
A itself while one extremum point appears for the first 
derivative of A. In view of the small 4x4 cluster used in 
our calculation, the good agreement of our results with 
previous studies is truly impressive. In other words, we 
demonstrate that the polarization spectroscopy method 
can be used to determine the quantum phase transition 
points straightforwardly and characterize distinct quan- 
tum phases. 

For the checkerboard antiferromagnet model, we cal- 
culate the noise signal A and its first derivative with 
respect to J 2 / ' J\. As shown in Fig. 4, there appears no 
extreme point for A itself, the evolution of A as a func- 
tion of J2/J1 can be roughly classified into three differ- 
ent regions, two flatter regions for both J2/J1 < 0.8 and 
J2/J1 > 1.2, and one steeper region for 0.8 < J2/J1 < 
1.2. By analyzing the first derivative A as a function of 
J2I J\, there is a clear maximum point around J2/J1 = 1. 
From the previous studies [2i| Ho] , there might be three 



different phases dividing by two quantum phase transi- 
tions points which are located at {J2/ Ji)c ~ 0.8 and 1.1, 
respectively. In this case, there is no good one-to-one 
correspondence between the extreme points of noise sig- 
nal and the quantum phase transition points. Since the 
unit cell of checkerboard lattice is larger than that of 
simple square lattice, a larger cluster is required to do 
numerical calculations. Unfortunately, it is now beyond 
our computational capabilities. 



V. SUMMARY 

Simulation of geometrically frustrated quantum anti- 
ferromagnets of ultracold fermionic atoms in optical lat- 
tices and detection of many-body states become a very 
hot topic both experimentally and theoretically. In this 
paper, we have shown such strongly correlated systems 
can be simulated by a specific optical engineering tech- 
nique. In particular, both the frustrated J1-J2 Heisen- 
berg model and checkerboard antiferromagnet model in 
the strongly frustrated regime can be realized for ultra- 
cold fermionic atoms in optical lattices. Moreover we 
employ the polarization spectroscopy method to probe 
the quantum many-body states of such systems. This 
method is designed to detect the spatially varying spin- 
spin correlation functions. In particular, for the J1-J2 
model, the quantum phase transition points can be de- 
termined straightforwardly and distinct quantum phases 
can be well characterized. The polarization spectroscopy 
method is quite promising to probe various magnetic or- 
ders of geometrically frustrated systems. 
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